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ABSTRACT 

The  Princeton  Ocean  Model  (POM)  has  been  implemented  in  the  South  China  Sea  for  hindcast  of  circulation 
and  thermohaline  structure.  A  two-step  technique  is  used  to  initialize  POM  with  temperature,  salinity,  and  velocity 
for  1  April  1998  and  integrate  it  from  1  April  1998  with  synoptic  surface  forcing  for  3  months  with  and  without 
data  assimilation.  Hydrographic  and  current  data  acquired  from  the  South  China  Sea  Monsoon  Experiment 
(SCSMEX)  from  April  through  June  1998  are  used  to  verify,  and  to  assimilate  into,  POM.  The  mean  SCSMEX 
data  (Apr-Jun  1998)  are  about  0.5°C  warmer  than  the  mean  climatological  data  above  the  50-m  depth,  and 
slightly  cooler  than  the  mean  climatological  data  below  the  50-m  depth,  and  are  fresher  than  the  climatological 
data  at  all  depths  and  with  the  maximum  bias  (0.2-0.25  ppt)  at  75-m  depth. 

POM  without  data  assimilation  has  the  capability  to  predict  the  circulation  pattern  and  the  temperature  field 
reasonably  well,  but  has  no  capability  to  predict  the  salinity  field.  The  model  errors  have  Gaussian-type  distri¬ 
bution  for  temperature  hindcast,  and  non-Gaussian  distribution  for  salinity  hindcast  with  six  to  eight  times  more 
frequencies  of  occurrence  on  the  negative  side  than  on  the  positive  side.  Data  assimilation  enhances  the  model 
capability  for  ocean  hindcast,  if  even  only  conductivity-temperature-depth  (CTD)  data  are  assimilated.  When 
the  model  is  reinitialized  using  the  assimilated  data  at  the  end  of  a  month  (30  Apr;  31  May  1998)  and  the  model 
is  run  for  a  month  without  data  assimilation  (hindcast  capability  test),  the  model  errors  for  both  temperature 
and  salinity  hindcast  are  greatly  reduced,  and  they  have  Gaussian-type  distributions  for  both  temperature  and 
salinity  hindcast.  Hence,  POM  gains  capability  in  salinity  hindcast  when  CTD  data  are  assimilated. 


1.  Introduction 

The  South  China  Sea  (SCS)  is  a  semienclosed  tropical 
sea  located  between  the  Asian  landmass  to  the  north 
and  west,  the  Philippine  Islands  to  the  east,  Borneo  to 
the  southeast,  and  Indonesia  to  the  south  (Fig.  1),  a  total 
area  of  3.5  X  106  km2.  It  includes  the  shallow  Gulf  of 
Thailand  and  connections  to  the  East  China  Sea  (through 
the  Taiwan  Strait),  the  Pacific  Ocean  (through  the  Luzon 
Strait),  the  Sulu  Sea,  the  Java  Sea  (through  the  Gasper 
and  Karimata  Straits),  and  to  the  Indian  Ocean  (through 
the  Strait  of  Malacca).  All  of  these  straits  are  shallow 
except  Luzon  Strait  whose  maximum  depth  is  1800  m. 
Consequently,  the  SCS  is  considered  a  semienclosed 
sea.  The  complex  topography  includes  a  broad  shallow 
shelf  in  the  south-southwest;  the  continental  shelf  of  the 
Asian  landmass  in  the  north,  extending  from  the  Gulf 
of  Tonkin  to  Taiwan  Strait;  a  deep,  elliptical  shaped 
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basin  in  the  center;  and  numerous  reef  islands  and  un¬ 
derwater  plateaus  scattered  throughout.  The  shelf  that 
extends  from  the  Gulf  of  Tonkin  to  the  Taiwan  Strait  is 
consistently  nearly  70  m  deep,  and  averages  150  km  in 
width;  the  central  deep  basin  is  1900  km  along  its  major 
axis  (northeast-southwest)  and  approximately  1100  km 
along  its  minor  axis,  and  extends  to  over  4000  m  deep. 
The  south-southwest  SCS  shelf  is  the  submerged  con¬ 
nection  between  Southeast  Asia,  Malaysia,  Sumatra, 
Java,  and  Borneo  and  reaches  100-m  depth  in  the  mid¬ 
dle;  the  center  of  the  Gulf  of  Thailand  is  about  70  m 
deep. 

The  SCS  is  subjected  to  a  seasonal  monsoon  system 
(Wyrtki  1961).  From  April  to  August,  the  weaker  south¬ 
westerly  summer  monsoon  winds  result  in  a  monthly 
mean  wind  stress  of  just  over  0.1  N  m~2  (Fig.  2a).  From 
November  to  March,  the  stronger  northeasterly  winter 
monsoon  winds  corresponds  to  a  maximum  monthly 
mean  wind  stress  of  nearly  0.3  N  m2  (Fig.  2b).  The 
transitional  periods  are  marked  by  highly  variable  winds 
and  surface  currents. 

A  survey  by  Wyrtki  (1961),  which  was  modeled  by 
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Fig.  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the  South  China  Sea.  Numbers  show 

the  depth  (m). 


Chu  et  al.  (1999a),  using  the  Princeton  Ocean  Model 
(POM)  revealed  complex  temporal  and  spatial  features 
of  the  surface  currents  in  both  the  SCS  and  the  sur¬ 
rounding  waters.  The  surface  circulations  for  the  sum¬ 
mer  (Jun-Aug)  show  the  strong  western  boundary  cur¬ 
rent  flowing  northward  to  northeastward  (Fig.  3a).  By 
April,  the  northeast  monsoon  has  decayed  and  the  main 
SCS  current  consisted  of  two  large  areas  of  anticyclonic 
flow  in  the  western  SCS,  one  in  the  north  and  the  other 
in  the  south.  With  the  beginning  of  the  southwest  mon¬ 
soon  in  May  and  June,  a  wind  drift  current  is  formed. 
Off  the  coast  of  Vietnam,  westward  intensification  of 


this  current  is  clearly  visible.  The  wide,  uniform  drift 
in  the  northern  SCS  shows  a  deflection  of  the  current 
to  the  right  of  the  wind  (Fig.  3a).  The  winter  circulations 
indicate  a  reverse  pattern:  the  western  boundary  current 
flowing  southward  to  southeastward.  The  southern  SCS 
is  occupied  by  a  cyclonic  flow  (Fig.  3b).  The  simulated 
SCS  surface  circulation  patterns  agree  with  the  obser¬ 
vational  results  quite  well  (see  Plates  3  and  6  in  Wyrtki 
1961). 

The  observed  circulation  patterns  of  the  intermediate 
to  upper  layers  of  the  SCS  are  primarily  forced  by  the 
local  monsoon  systems  (Wyrtki  1961),  with  contribu- 
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Fig.  2.  Climatological  wind  stress  vectors  obtained  from  COADS  (da  Silva  1994):  (a) 

Jun  and  (b)  Dec. 


tions  from  the  Kuroshio  Current  via  the  Bashi  Channel. 
The  Kuroshio  entered  the  SCS  through  the  southern  side 
of  the  channel,  then  executed  a  tight,  anticyclonic  turn, 
and  exited  the  SCS  near  Taiwan.  An  estimated  6.5  Sv 
(1  Sv  =  106  m3  s_1)  of  the  intrusion  passed  through 
Bashi  Channel  (Chu  and  Li  2000)  using  the  navy’s  cli¬ 
matological  temperature  and  salinity  dataset  (public  do¬ 
main)  with  1/2°  X  1/2°  resolution  by  the  P-vector  meth¬ 
od.  This  flow  exerted  a  strong  influence  on  the  properties 
of  the  northern  SCS  waters  and  is  believed  to  contribute 
to  currents  in  the  Taiwan  Strait  (Hu  and  Liu  1992). 


The  Naval  Oceanographic  Office  (NAVOCEANO) 
and  the  Fleet  Numerical  Meteorology  and  Oceanogra¬ 
phy  Center  (FNMOC)  are  the  major  operational  envi¬ 
ronmental  processing  centers.  Both  NAVOCEANO  and 
FNMOC  are  interested  in  comparisons  of  advanced  nu¬ 
merical  analysis  and  ocean  nowcast/hindcast  models 
leading  to  model  improvements.  The  South  China  Sea 
Monsoon  Experiment  (SCSMEX)  provides  a  unique  op¬ 
portunity  for  such  an  evaluation.  Hydrographic  and  cur¬ 
rent  data  acquired  from  SCSMEX  for  April  through  June 
1998  are  used  to  verify,  and  to  assimilate  into,  POM. 


(a) 


(b) 


Fig.  3.  Mean  surface  circulation  simulated  using  POM-SCS:  (a)  summer  and  (b) 
winter.  [After  Chu  et  al.  (1999a).] 
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SCS  CTD  Data  Station:  1742 


SCS  ADCP  Data  Station:  33923 
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Fig.  4.  The  SCSMEX  oceanographic  observations:  (a)  CTD  stations  and  (b)  ADCP  stations. 


2.  Oceanographic  observations  during  SCSMEX 

a.  Description 

SCSMEX  is  a  large-scale  experiment  to  study  the 
water  and  energy  cycles  of  the  Asian  monsoon  regions 
with  the  goal  (SCSMEX  Science  Working  Group  1995) 
to  provide  a  better  understanding  of  the  key  physical 
processes  responsible  for  the  onset,  maintenance,  and 
variability  of  the  monsoon  over  Southeast  Asia  and 
southern  China  leading  to  improved  predictions.  It  in¬ 
volves  the  participation  of  all  major  countries  and  re¬ 
gions  of  East  and  Southeast  Asia,  as  well  as  Australia 
and  the  United  States. 

SCSMEX  had  both  atmospheric  and  oceanic  com¬ 
ponents.  The  oceanic  intensive  observational  period 
(IOP)  was  from  April  through  June  1998  with  shipboard 
measurements,  ATLAS  moored  array,  and  drifters.  The 
hydrographic  data  collected  during  the  SCSMEX  IOP 
went  through  quality  control  procedures,  such  as  a  min¬ 
imum-maximum  check  (e.g.,  disregarding  any  temper¬ 
ature  data  less  than  — 2°C  and  greater  than  40°C),  an 
error  anomaly  check  (e.g.,  rejecting  temperature  data 
deviating  more  than  7°C  from  climatology),  a  ship¬ 
tracking  algorithm  (screening  out  data  with  obvious  ship 
position  errors),  a  maximum  number  limit  (limiting  a 
maximum  number  of  observations  within  a  specified  and 
rarely  exceeded  space-time  window),  and  a  buddy 
check  (tossing  out  contradicting  data).  The  climatolog¬ 


ical  data  used  for  quality  control  are  depicted  in  Chu  et 
al.  (1997a,b).  After  the  quality  control,  the  SCSMEX 
oceanographic  dataset  contains  1742  CTD  stations  (Fig. 
4a)  and  33  923  Acoustic  Doppler  Current  Profiler 
(ADCP)  stations  (Fig.  4b).  The  majority  of  the  CTDs 
were  nominally  capable  of  reaching  a  maximum  depth 
of  2000  m. 

b.  T-S  diagram 

The  monthly  CTD  station  distribution  during 
SCSMEX  1998  (Fig.  5)  shows  good  coverage  for  the 
SCS  basin.  The  corresponding  T-S  diagrams  (opposite- 
S  shape  T-S  curves)  clearly  show  the  existence  of  five 
water  masses:  the  SCS  surface  water  (SCSSW,  warm 
and  less  fresh),  the  SCS  subsurface  water  (SCSSSW, 
less  warm  and  salty),  the  SCS  intermediate  water 
(SCSIW,  less  cool  and  fresh),  the  SCS  deep  water 
(SCSDW,  cool  and  fresher),  and  the  SCS  coastal  water 
(SCSCW,  cool  and  fresh).  The  characteristics  of  the  five 
water  masses  are  illustrated  in  Table  1. 

c.  Current  field 

The  SCSMEX  ADCP  data  at  20-m  depth  show  the 
existence  of  multieddy  structure  in  SCS  (Fig.  6a)  with 
north  cyclonic-south  anticyclonic  dipoles  (9°-14°N, 
1 10°-1 13°E)  near  the  Vietnamese  Bight,  an  anticyclonic 
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Fig.  5.  CTD  stations  during  SCSMEX-1998  and  corresponding  T-S  diagrams. 
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Table  1.  Hydrographic  features  of  the  five  SCS  waters. 


Water  mass 

Temperature  (°C) 

Salinity  (ppt) 

Depth  (m) 

scssw 

27.5 

33.75 

<50 

scsssw 

20.5 

34.72 

50-200 

SCSIW 

7.0 

34.36 

200-600 

SCSDW 

2.0 

34.64 

>1000 

scscw 

17.0 

33.24 

<50 

Table  2.  Bimonthly  variation  of  volume  transport  (Sv)  at  the  lateral 
open  boundaries.  The  positive  (negative)  values  mean  outflow  (in¬ 
flow)  and  were  taken  from  Wyrtki  (1961). 


Feb 

Apr 

Jun 

Aug 

Oct 

Dec 

Gaspar  and  Kari- 
mata  Straits 

4.4 

0.0 

-4.0 

-3.0 

1.0 

4.3 

Luzon  Strait 

-3.5 

0.0 

3.0 

2.5 

-0.6 

-3.4 

Taiwan  Strait 

-0.9 

0.0 

1.0 

0.5 

-0.4 

-0.9 

eddy  (8°-13°N,  115°-120°E)  west  of  Palawan  Island, 
and  a  cyclonic  eddy  (17°-21°N,  110°-113°E)  east  of 
Hainan  Island. 


3.  Model 

a.  Description 


where  KM  is  the  vertical  mixing  coefficient  for  momen¬ 
tum,  and  (u,  v)  and  (t0jc,  t0  )  are  the  two  components 
of  the  water  velocity  and  wind  stress  vectors,  respec¬ 
tively. 

Surface  thermal  forcing  is  determined  by 


POM  is  a  time-dependent,  primitive  equation  circu¬ 
lation  model  on  a  three-dimensional  grid  that  includes 
realistic  topography  and  a  free  surface.  Developed  at 
Princeton  University  (Blumberg  and  Mellor  1987),  the 
model  was  specifically  designed  to  accommodate  me- 
soscale  phenomena,  including  the  often  nonlinear  pro¬ 
cesses  (such  as  upwelling  and  eddy  dynamics)  com¬ 
monly  found  in  estuarine  and  coastal  oceanography. 
Several  years  ago  POM  was  implemented  in  the  domain 
(3.06°S-25.07°N,  98.84°-121.16°E)  covering  the  whole 
SCS  and  surrounding  land  and  islands  (Chu  et  al.  1998a, 
1999a).  POM  has  the  capability  to  simulate  the  SCS 
multieddy  structure  detected  from  both  historical  (Chu 
et  al.  1997a,b)  and  synoptic  data  (Chu  et  al.  1998c)  and 
to  simulate  the  SCS  response  to  a  tropical  cyclone  (Chu 
et  al.  2000b).  POM  is  also  a  useful  tool  for  investigating 
the  mechanisms  for  the  seasonal  variability  (Chu  et  al. 
1999a)  and  eddy  formation  (Chu  et  al.  1998a).  Recently, 
a  coastal  atmosphere-ocean  coupled  system  was  de¬ 
veloped  with  POM  as  the  oceanic  component  and  the 
Pennsylvania  State  University-National  Center  for  At¬ 
mospheric  Research  fifth-generation  Mesoscale  as  the 
atmospheric  component  (Chu  et  al.  1999b). 

In  this  study,  the  horizontal  spacing  is  0.179°  by 
0.175°  (approximately  20-km  resolution)  and  there  are 
23  vertical  sigma  coordinate  levels.  Consequently,  the 
model  contains  125  X  162  X  23  fixed  grid  points.  The 
model  uses  realistic  bathymetric  data  from  the  NA- 
VOCEANO  Digital  Bathymetry  Data  Base  (NAVO- 
CEANO  1998)  5  min  by  5  min  resolution.  The  hori¬ 
zontal  diffusion  coefficient  is  chosen  to  be  0.2  for  this 
application.  The  bottom  stress  r(,  is  assumed  to  follow 
a  quadratic  law  (Blumberg  and  Mellor  1987). 


b.  Atmospheric  forcing 

The  atmospheric  forcing  for  the  SCS  application  of 
POM  includes  mechanical  and  thermohaline  forcing. 
The  wind  forcing  is  determined  by 


PoKm\ 


du  dv 
dz’  dz 


=  (Tte,  Toy), 


(1) 


v  —  =  (  l  4-  rtn 

Kn  ^  (X~^OBS 


e) 


(2) 


^7  =  «I  Qs  +  ot2C(SOBS  -  S), 
dz 


Qs  =  S(E  -  Pr),  (3) 

where  0OBS  and  SOBS  are  the  observed  potential  temper¬ 
ature  and  salinity,  KH  and  Ks  are  the  vertical  mixing 
coefficients  for  temperature  and  salinity,  C  is  the  spe¬ 
cific  heat,  Qh  and  Qs  are  surface  net  heat  and  salinity 
fluxes,  and  E  and  P,  are  evaporation  and  precipitation 
rates,  respectively.  The  relaxation  coefficient  C  is  the 
reciprocal  of  the  restoring  time  period  for  a  unit  volume 
of  water.  The  two  weights,  a,  =  1  and  a2  =  0,  would 
specify  only  flux  forcing  is  applied;  cr,  =  0  and  a2  = 
1  would  specify  that  only  restoring-type  forcing  is  ap¬ 
plied.  The  relaxation  coefficient  C  is  taken  to  be  0.7  m 
day-1,  which  is  equivalent  to  a  relaxation  time  of  43 
days  for  an  upper  layer  30  m  thick  (Chu  et  al.  1998a, 
1999a-c).  The  net  effect  is  to  prevent  any  large  devi¬ 
ations  from  climatology.  Chu  et  al.  (1998b)  pointed  out 
that  the  pure  restoring  boundary  condition  (a,  =  0  and 
a2  =  1)  is  not  applicable  to  coastal  regions. 


c.  Lateral  boundary  forcing 

Rigid  lateral  boundaries,  that  is,  the  modeled  ocean 
bordered  by  land,  were  defined  using  a  free-slip  con¬ 
dition  for  velocity  and  a  zero  gradient  condition  for 
temperature  and  salinity.  No  advective  or  diffusive  heat, 
salt,  or  velocity  fluxes  occur  through  these  boundaries. 
Open  boundaries,  where  the  numerical  grid  ends  but  the 
fluid  motion  is  unrestricted,  were  treated  as  radiative 
boundaries.  When  the  water  flows  into  the  model  do¬ 
main,  temperature  and  salinity  at  the  open  boundary  are 
likewise  prescribed  from  the  climatological  data  (Lev- 
itus  1982).  When  water  flows  out  of  the  domain,  the 
radiation  condition  is  applied, 

£(0,  5)  +  t/„f  (0,  S)  =  0, 

dt  dn 


z=0 


(4) 
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Fig.  6.  Horizontal  velocity  vectors  (May-Jun  1998)  at  20-m  depth  (a)  obtained  from  ADCP  observations  and  (b)  predicted  from 

MD2. 


where  the  subscript  n  is  the  direction  normal  to  the 
boundary. 

The  volume  transports  through  major  straits  of  the 
SCS  are  uncertain.  For  example,  different  authors  re¬ 
ported  volume  transport  through  the  Luzon  Strait  rang¬ 
ing  from  8  to  10  Sv  (e.g.,  Huang  et  al.  1994),  6.5  Sv 
(Chu  and  Li  2000)  to  2-3  Sv  (Wyrtki  1961;  Metzger 
and  Hurlburt  1996).  It  is  not  clear  how  to  choose  one 
among  various  estimates  of  the  volume  transport.  De¬ 
spite  being  old,  Wyrtki’s  (1961)  data  provide  a  balanced 
estimate  of  volume  transports  for  the  Luzon  Strait,  Tai¬ 
wan  Strait,  and  Gasper  and  Karimata  Straits  with  sea¬ 
sonal  variations  (Table  2).  Since  there  are  no  reliable 
estimates  at  Balabac  Channel,  Mindoro  Strait,  and  the 
Strait  of  Malacca,  we  assumed  zero  transport  there.  Such 
a  treatment,  especially  at  Mindoro  Strait,  may  distort 
the  solution. 

d.  Two-step  initialization 

First,  we  integrate  POM  for  3  yr  from  zero  velocity 
and  January  climatological  temperature  and  salinity 
fields  (Levitus  1982)  with  climatological  monthly  mean 
surface  wind  stress  from  the  Comprehensive  Ocean  and 
Atmosphere  Data  Set  (COADS;  da  Silva  et  al.  1994) 
and  restoring-type  thermohaline  surface  forcing  (a,  = 
0,  a2  =  1),  which  are  relaxed  to  the  surface  monthly 
mean  values.  Second,  we  continue  to  integrate  POM  for 


2  yr  and  3  months  with  monthly  mean  surface  wind 
stress  and  heat  and  salt  fluxes  (al  =  1 ,  a2  =  0)  from 
COADS.  The  final  states  of  temperature,  salinity,  and 
velocity  fields  are  taken  as  the  initial  ocean  conditions 
(1  Apr  1998)  for  model  evaluation. 

e.  CTD  data  assimilation 

We  use  the  intermittent  optimal  interpolation  (OI) 
scheme  developed  at  Harvard  University  (Robinson  et 
al.  1996;  Lozano  et  al.  1996)  to  assimilate  the  SCSMEX 
CTD  data  daily.  The  large-scale  OI  was  used  to  estimate 
the  background  mean  with  a  decorrelation  scale  of  450 
km.  The  mesoscale  OI  was  used  to  map  the  observa¬ 
tional  anomaly  from  the  background  mean  field  into  a 
regular  grid  with  a  spatial  decorrelation  scale  of  75  km 
and  a  temporal  decorrelation  scale  of  10  days  (Chu  et 
al.  1998c).  The  spatial  and  temporal  decorrelation  scales 
associated  with  the  mean  field  were  estimated  from  the 
covariance  matrix  obtained  from  the  navy’s  Master 
Oceanographic  Observational  Data  Set  for  the  SCS  (see 
Chu  et  al.  1997c).  The  large-scale  mean  was  computed 
using  the  temperature  and  salinity  profiles.  The  decor¬ 
relation  length  scale  for  the  fluctuations  from  the  mean 
was  estimated  from  the  synoptic  data  (Chu  et  al.  1998c). 
The  decorrelation  timescale  was  selected  to  ensure  syn- 
opticity.  The  SST  horizontal  scales  seem  somewhat  larg¬ 
er  than  those  of  temperatures  in  the  seasonal  thermo- 
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Fig.  7.  Scatter  diagram  of  the  MD,  data  the  SCSMEX  data. 
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Fig.  9.  The  rmse  between  the  MD,  and  the  SCSMEX  CTD  data  (temperature,  °C;  salinity,  ppt). 
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Fig.  10.  The  bias  between  the  MD,  and  the  SCSMEX  CTD  data  (temperature,  °C;  salinity,  ppt). 
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Fig.  11.  Skill  score  for  POM  without  data  assimilation  (MD,  data). 
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Fig.  12.  Scatter  diagram  of  the  MD2  data  the  SCSMEX  data. 


cline.  In  this  study,  we  choose  horizontal  scales  to  be 
uniform  in  the  vertical  for  convenience. 

/.  Model  integration 

We  integrate  the  POM-SCS  model  from  the  state  of 
1  April  1998  with  daily  surface  wind  stress  and  heat 
and  salt  fluxes  (a,  =  1 ,  a2  =  0)  for  3  months  from 
National  Centers  for  Environmental  Prediction  data  in 
two  cases:  1)  without  any  data  assimilation  (defined  as 
MD,)  and  2)  with  daily  SCSMEX-CTD  data  assimi¬ 
lation  (defined  as  MD,). 

To  verify  the  model’s  hindcast  capability,  we  should 
define  hindcast  periods  during  which  the  observational 
data  are  not  assimilated.  The  hindcast  period  for  MD, 
is  3  months  (1  Apr-30  Jun  1998);  MD,  was  evaluated 
simply  by  direct  comparison  between  the  modeled  and 
SCSMEX  data. 

The  hindcast  periods  for  MD,  are  May  and  lune  1998. 


For  the  May  (Jun)  MD,  hindcast,  we  use  the  states  of 
30  April  (31  May)  from  the  assimilation  run  (MD,)  as 
the  new  initial  conditions  to  integrate  POM  for  a  month 
with  no  data  assimilation.  Then  MD2  was  evaluated  by 
comparison  between  the  reinitialization  model  and  the 
SCSMEX  CTD  data  for  May  and  June  1998.  If  MD, 
has  better  capability  than  MD,,  we  may  confirm  the 
value  added  by  the  CTD  data  assimilation. 

The  circulation  and  the  thermohaline  structures  from 
both  runs  (MD,  and  MD2)  show  similar  patterns.  The 
predicted  SCS  surface  circulation  is  generally  anticy- 
clonic  during  the  SCSMEX  period  (Apr-Jun  1998),  with 
a  strong  western  boundary  current,  splitting  and  par¬ 
tially  leaving  the  coast;  the  bifurcation  point  is  at  14°N 
in  May,  and  shifts  southward  in  June.  Recently,  we  have 
reported  POM-SCS  model  results  (e.g.,  Fig.  3)  in  a 
series  of  papers  (Chu  et  al.  1999a-c;  Chu  et  al.  2000a, b), 
and  thus  in  this  paper,  we  mainly  show  the  results  of 
the  model  verification. 
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Fig.  13.  Histogram  of  model  (MD2)  errors  of  temperature  (°C)  and  salinity  (ppt). 


The  predicted  20-m  depth  (Fig.  6b)  velocity  vector 
fields  (averaged  over  May-Jun  1998)  from  MD,  show 
the  multieddy  structure  with  a  central  SCS  anticyclonic 
gyre  surrounded  by  north  cyclonic-south  anticyclonic 
dipoles  (9°-14°N,  110°-113°E)  near  the  Vietnamese 
Bight,  an  anticyclonic  eddy  (8°-13°N,  1 15°-120°E)  west 
of  Palawan  Island,  and  a  cyclonic  eddy  (17°-21°N, 
110°-113°E)  east  of  Hainan  Island.  The  predicted  mul¬ 
tieddy  structure  is  consistent  with  the  SCSMEX  ADCP 
data. 

4.  Methodology  of  verification 

The  observational  data  are  located  at  depth  z.  We 
interpolate  the  model  data  into  the  observational  points 
(xh  ypz,  t )  and  form  two  sets  of  modeled  data:  MD,  and 
MD,. 


a.  Error  distribution 

The  difference  of  the  modeled  and  observed  vari¬ 
able  i/f, 

Ai//(v,  yP  z,  t)  =  I yp  z,  t) 

-  i li0{x„  yP  z,  t ),  (5) 

represents  the  model  error.  Here  t j/m  and  i //„  are  the 
variables  from  model  and  observation,  respectively. 
We  may  take  the  probability  histogram  of  A  iff  as  the 
error  distribution. 

b.  Basic  statistical  parameters 

Three  parameters:  Bias,  mean-square  error  (mse), 
and  root-mean-square  error  (rmse). 
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Fig.  14.  The  rmse  between  the  MD,  and  the  SCSMEX  CTD  data  (temperature-  °C;  salinity,  ppt). 


biasfc  o  =  ^  E  E  At//(x,.,  yp  z,  t), 

Jy  i  j 

msefc  0  =  EE  ^[Ai//(-T,  yjt  z,  OP, 


rmse(z,  t )  =  Vmse(z,  t),  (6) 

are  commonly  used  for  evaluation  of  the  model  per¬ 
formance.  Here  N  is  the  total  number  of  horizontal 
points,  and  [i|i,„(z,  t),  i|i„(z,  0]  and  [c rm(z,  0,  ct0(z,  0] 
are  horizontal  means  and  standard  deviations  of  the 
modeled  and  observed  data. 


c.  Model  accuracy  and  skill 

Model  accuracy  is  usually  defined  as  the  average 
degree  of  correspondence  between  modeled  and  ob¬ 


servational  data.  Thus,  the  mse  or  rmse  represents  a 
prototype  measure  of  accuracy.  Model  skill,  on  the 
other  hand,  is  defined  as  the  model  accuracy  relative 
to  the  accuracy  of  hindcasts  produced  by  some  ref¬ 
erence  procedure  such  as  climatology  or  persistence. 
To  measure  the  model  skill,  we  may  compute  the  re¬ 
duction  of  mse  over  the  climatological  hindcasts 
(Murphy  1988), 


SS  =  1  - 


mse(m,  o) 
mse(c,  o)  ’ 


(7) 


where  mse  (m,  o)  and  mse(c,  o )  are  the  mse  of  the 
model  and  climatological  hindcasts,  respectively.  The 
parameter  SS  is  called  the  skill  score,  which  is  pos¬ 
itive  (negative)  when  the  accuracy  of  the  hindcasts  is 
greater  (less)  than  the  accuracy  of  the  reference  hind¬ 
casts  (climatology).  Moreover,  SS  =  1  when  mse(m. 
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o)  =  0  (perfect  hindcasts)  and  SS  =  0  when  mse(m, 
o)  =  mse(c,  o).  To  compute  mse(c,  o),  we  interpolate 
the  climatological  monthly  temperature  and  salinity 
data  (Levitus  1982)  into  the  observational  points  (x,, 
yj,  z,  t ). 

5.  Evaluation  of  POM  without  data  assimilation 

We  compare  the  MD,  data  against  the  SCSMEX 
CTD  data  to  verify  POM-SCS  model  without  data 
assimilation. 

a.  Statistical  evaluation 

The  easiest  way  to  verify  POM-SCS  performance 
is  to  plot  the  MD,  data  against  SCSMEX  CTD  data 
(Fig.  7).  The  scatter  diagrams  for  temperature  show 
the  points  clustering  around  the  line  of  Tm  =  Ta.  The 
scatter  diagrams  for  salinity  show  much  more  spread¬ 
ing  out  of  the  points  around  the  line  Sm  =  S„.  This 
result  indicates  better  performance  in  temperature 
hindcast  than  in  salinity  hindcast. 

The  model  errors  for  temperature  hindcast  have 
Gaussian-type  distribution  with  mean  values  approx¬ 
imately  zero  and  values  of  standard  deviation  (STD) 
increasing  with  integration  time  from  1 .48°C  in  April 
and  1.63°C  in  June  (Fig.  8),  but  the  model  errors  for 
salinity  hindcast  reveal  non-Gaussian  distribution 
(Fig.  8)  with  quite  steady  STD  values  (0.186-0.189 
ppt)  and  with  six  to  eight  times  more  frequencies  of 
occurrence  on  the  negative  side  than  on  the  positive 
side:  3318  (-)  versus  576  (  +  )  in  April,  20  527  (-) 
versus  3224  (  +  )  in  May,  and  74  707  (— )  versus  1672 
(+)  in  June.  These  results  indicate  that  the  model 
underpredicts  the  salinity. 

b.  Error  estimation 

The  rmse  of  temperature  (Fig.  9)  between  the  MD, 
data  and  the  SCSMEX  data  in  April  is  around  1°C 
from  the  surface  to  350  m  and  decreases  with  depth 
to  near  0.5°C  at  500-m  depth.  The  value  of  rmse(m, 
o )  increases  with  time.  For  the  June  hindcast,  rmse(m, 
o)  is  around  1.3°C  in  the  upper  layer  from  the  surface 
to  250  m,  and  the  value  of  rmse(m,  o )  is  greater  than 
the  value  of  rmse(c,  o)  at  most  of  the  depths.  This 
result  indicates  the  POM-SCS  model  loses  its  pre¬ 
dictability  for  temperature  in  2  months  from  the  initial 
state  (1  Apr  1998)  without  data  assimilation. 

The  rmse  of  salinity  (Fig.  9)  between  the  MD,  data 
and  the  SCSMEX  data  is  smaller  than  that  between 
the  climatological  data  and  the  SCSMEX  data  only 
in  the  shallow  upper  layer  (surface  to  60-m  depth)  in 
April,  and  greater  than  that  between  the  climatolog¬ 
ical  data  and  the  SCSMEX  data  elsewhere  in  April 
and  all  the  depths  in  May  and  June. 

The  mean  SCSMEX  data  in  April-June  1998  are 
around  0.5°  warmer  than  the  mean  climatological  data 


above  the  50-m  depth  and  colder  than  the  mean  cli¬ 
matological  data  below  the  50-m  depth,  and  are  fresh¬ 
er  than  the  climatological  data  in  all  depths,  with  the 
maximum  bias  (0.2-0.25  ppt)  at  75-m  depth  (Fig.  10). 

The  POM-SCS  modeled  mean  temperature  is  about 
0.5°C  (0.5°-1.0°C)  cooler  (warmer)  than  the 
SCSMEX  data  above  (below)  50-m  depth,  and  the 
POM-SCS  modeled  mean  salinity  is  less  that  the 
SCSMEX  data  with  the  maximum  fresh  bias  of  0.25 
ppt  in  depths  between  50  and  100  m  (Fig.  10). 

The  skill  score  of  the  temperature  hindcast  (Fig. 
11)  is  positive  in  April  and  May,  and  becomes  neg¬ 
ative  in  June,  which  indicates  that  the  model  loses  the 
capability  to  predict  the  temperature  in  2  months.  The 
skill  score  of  the  salinity  hindcast  (Fig.  1 1)  is  negative 
almost  everywhere,  which  indicates  the  model  has 
little  capability  to  predict  salinity. 

6.  Evaluation  of  POM  with  the  CTD  data 
assimilation 

We  compare  the  MD,  data  against  the  SCSMEX 
CTD  data  to  verify  the  valueadded  of  reinitialization 
using  data  assimilation. 

a.  Statistical  evaluation 

We  plot  the  MD,  hindcast  data  against  SCSMEX 
CTD  data  (Fig.  12).  The  points  in  both  temperature 
and  salinity  scatter  diagrams  are  more  clustered 
around  the  lines  Tm  =  Ta  and  Sm  =  S0  in  the  assim¬ 
ilated  model  (Fig.  12)  than  in  the  nonassimilated  mod¬ 
el  (Fig.  7).  This  result  indicates  the  improvement  of 
model  reinitialization  using  the  SCSMEX  CTD  data 
in  an  ocean  hindcast.  For  both  a  temperature  and  sa¬ 
linity  hindcast,  the  model  errors  have  Gaussian-type 
distributions  (Fig.  13)  with  much  smaller  STD  values 
(e.g.,  1. 19°C,  0. 143  ppt  in  Jun)  compared  to  MD,  (Fig. 
8).  However,  the  salinity  hindcast  errors  are  two  to 
three  times  greater  on  the  negative  side  than  on  the 
positive  side:  13  687  (— )  versus  7111  (  +  )  in  May, 
and  1 1  643  (  — )  versus  3414  (  +  )  in  June.  These  results 
indicate  that  the  model  (MD,)  still  underpredicts  the 
salinity. 

b.  Error  estimation 

The  rmse  for  temperature  (Fig.  14)  is  0.2°C  at  the 
surface,  increases  with  depth  to  a  maximum  value  of 
1.2°C  (May  1998)  or  of  1.4  (Jun  1998)  at  50-100-m 
depth,  and  then  decreases  with  depth  to  a  minimum 
value  of  0.3°C  (May  1998)  or  0.5°C  at  500  m.  The 
mean  rmse  is  around  0.6°C.  The  rmse  for  salinity  (Fig. 
14)  is  near  zero  at  the  surface,  increases  with  depth 
to  a  maximum  value  of  0.22  ppt  (May  1998)  or  of 
0.21  ppt  (Jun  1998)  at  25-75-m  depth,  and  then  de¬ 
creases  with  depth  to  a  minimum  value  of  0.03  ppt 
at  500  m.  The  mean  rmse  is  around  0.06  ppt.  Both 
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Fig.  15.  The  bias  between  the  MD,  and  the  SCSMEX  CTD  data  (temperature,  °C;  salinity,  ppt). 


temperature  and  salinity  rmse’s  between  the  MD,  data 
and  the  SCSMEX  data  (Fig.  14)  are  smaller  than  that 
between  the  MD,  data  and  the  SCSMEX  data  (Fig. 
9),  especially  for  salinity.  The  rmse(m,  o)  for  salinity 
becomes  smaller  than  rmse(c,  o )  for  almost  all  the 
depths  in  May  and  June.  These  results  indicate  that 
POM-SCS  gains  the  capability  of  salinity  hindcast  if 
we  reinitialize  the  model  monthly  with  the  SCSMEX 
CTD  data.  The  values  of  rms e(m,  o )  for  temperature 
are  also  reduced,  compared  to  the  MDj  case. 

The  POM-SCS  (with  data  assimilation)  modeled 
mean  temperature  is  about  0.2°C  warmer  in  May  (and 
about  0.5°C  warmer  in  Jun)  than  the  SCSMEX  data 
at  100-m  depth,  and  the  POM-SCS  (with  data  assim¬ 
ilation)  modeled  mean  salinity  is  less  than  the 
SCSMEX  data  with  the  maximum  fresh  bias  of  0.17 


ppt  at  20-m  depth  in  May  1998  and  0.1  ppt  in  depths 
between  25  and  75  m  (Fig.  15). 

The  skill  scores  for  both  temperature  and  salinity 
hindcasts  (Fig.  16)  are  much  improved  compared  to 
MD,  (Fig.  11),  especially  the  salinity  hindcast.  The 
skill  score  of  the  MD,  model  (Fig.  16)  for  the  salinity 
hindcast  is  positive  except  in  the  sublayer  (125-175 
m). 

7.  Conclusions 

1)  The  mean  SCSMEX  data  in  April-June  1998  are 
about  0.5°  warmer  than  the  mean  climatological 
data  (Fevitus  1982)  above  50-m  depth  and  cooler 
than  the  mean  climatological  data  below  50-m 
depth.  They  are  fresher  than  the  climatological 


1538 


Volume  18 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECHNOLOGY 


Temperature  (°C)  Salinity  (psu) 


Fig.  16.  Skill  score  for  POM  with  reinitialization  (MD,  data). 


data  at  all  depths,  with  the  maximum  bias  (0.2- 
0.25  ppt)  at  75-m  depth. 

2)  POM-SCS  with  pure  surface  flux  forcing  and 
without  data  assimilation  has  the  capability  to  pre¬ 
dict  the  temperature  filed  reasonably  well,  but  has 
little  capability  to  predict  the  salinity  field.  The 
model  errors  for  temperature  hindcast  have  a 
Gaussian-type  distribution,  but  the  model  errors 
for  salinity  hindcast  reveal  a  non-Gaussian  distri¬ 
bution  with  six  to  eight  time  more  frequencies  of 
occurrence  on  the  negative  side  than  on  the  pos¬ 
itive  side.  This  result  indicates  that  the  model  un¬ 
derpredicts  the  salinity.  The  skill  score  of  the  tem¬ 
perature  hindcast  is  positive  in  April  and  May,  and 
becomes  negative  in  June,  which  indicates  that  the 
model  will  lose  the  capability  to  predict  the  tem¬ 
perature  in  2  months  without  data  assimilation. 
The  skill  score  of  the  salinity  hindcast  is  negative 


almost  everywhere,  which  indicates  a  lack  of  ca¬ 
pability  in  salinity  hindcast. 

3)  Data  assimilation  enhances  the  model  capability 
for  an  ocean  hindcast,  even  if  only  CTD  data  are 
assimilated.  When  we  reinitialize  the  model  using 
the  assimilated  data  at  the  end  of  a  month  (30  Apr, 
30  May  1998),  the  model  errors  for  both  the  tem¬ 
perature  and  salinity  hindcasts  have  Gaussian-type 
distributions.  Both  temperature  and  salinity  rmse’s 
are  smaller  with  the  CTD  data  assimilation  than 
that  without  the  data  assimilation,  especially  for 
salinity.  With  the  data  assimilation,  the  mean  rmse 
for  temperature  is  around  0.6°C,  and  for  salinity 
is  around  0.06  ppt. 
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